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ABSTRACT 

(N| We investigate observational signatures of triggered star formation in bright rimmed 

clouds (BRCs) by using molecular line transfer calculations based on radiation- 
r ~T l hydrodynamic radiatively-driven-implosion models. We find that for BRCs the sepa- 

ration in velocity between the line profile peak of an optically thick and an optically 
Q^) thin line is determined by both the observer viewing angle and the density of the shell 

driving into the cloud. In agreement with observations, we find that most BRC line 
profiles are symmetric and that asymmetries can be either red or blue, in contrast to 
the blue-dominance expected for a collapsing cloud. Asymmetries in the line profiles 
arise when an optically thick line is dominated by the shell and an optically thin line is 
dominated by the cloud interior to the shell. The asymmetries are red or blue depend- 
ing on whether the shell is moving towards or away from the observer respectively. 
, ^ i Using the known motions of the molecular gas in our models we rule out the 'envelope 

expansion with core collapse' mechanism as the cause of the lack of blue-asymmetry in 
fN| our simulated observations. We show that the absence of a strong photon dominated 

region (PDR) around a BRC may not rule out the presence of triggered star formation: 
C^l if the BRC line profile has a strong blue component then the shell is expected to be 

driving towards the observer, suggesting that the cloud is being viewed from behind 
and the PDR is obstructed. This could explain why BRCs such as SFO 80, 81 and 86 
T "~J have a blue secondary peak and only a weak PDR inferred at 8 fim. Finally we also 

CO test the use of 12 CO, 13 CO and C 18 as diagnostics of cloud mass, temperature and 

column density. We find that the inferred conditions are in reasonable agreement with 
CO those from the models. Calculating the cloud mass assuming spherical symmetry is 

shown to introduce an error of an order of magnitude whereas integrating the column 
^ density over a given region is found to introduce an error of up to a factor of two. 
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1 INTRODUCTION 

Molecular line diagnostics are a widely used tool for investi- 
gating the conditions of astrophysical clouds and star forma- 
tion (e.g.|Schneps et al.|1980||De Vries et al. 12002] |Urquhart 



et al. 2006; N arayanan et al |2008||Rundle et al.|2010||Buckle" 



et al.[2012| Hatchell et al.[2013). Line profiles can yield in- 



formation about the kinematic motions of the molecular gas 



(e.g. 


Lee et al.|2004 


Tsamis et al.|2008 


Roberts et al.| 2010 


Stahler & Yen 2010 


Lou & Gao 2011 


) and ratios of line 



intensities can be used to infer the cloud properties such as 
optical depth, temperature, column density and mass (e.g. 
Myers et aL|l983| ). 



Bright rimmed clouds (BRCs) are objects which are be- 
lieved to be formed when shocks generated by nearby mas- 
sive stars drive into surrounding pre-existing density struc- 
tures, potentially triggering star formation in the radiatively 
driven implosion (RDI) scenario (e. g. |Sandford et al.||1982| 



Bertoldil[l989l JLefloch fe Lazareffl|1994| JKessel-Deynet fe| 
Burkert||2003| IGritschneder et al.| 2009) |Miao et al.| | 2009 | 



* E-mail: haworth@astro.ex.ac.uk 



Mackey Lim|2010[|Bisbas et al.|2011||Tremblin et al.|2012 
Haworth &; Harries||2012| and references therein). This is a 
highly kinematic process in which shock driving is expected 
to occur, potentially leading to the collapse of the cloud. 
As such, BRCs have been subjected to a large number of 
molecular line observations to try and identify RDI (e.g. 
Lefioch et aT]^997| |De Vries et aL] [20021 [Thompson et al.| 
2004| lUrquhart et al.|2006[|2009[|Morgan et al.||2009| ). 
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The conditions in the neutral gas of BRCs are frequently 
calculated based on the ratio of 13 CO to C 18 intensities 



following Myers et al. (1983). However alternative combi- 
nations of lines can also be used such as CS, HCO + , HCN 



and oth er CO isotopologues. [Lefloch et al. (1997) used CO, 
CS and|Thompson et al.|fl2004|) used x * 00,^*00 to calcu- 



late the conditions of BRCs in IC1848. They then compared 
the neutral cloud conditions with the ionized boundary layer 
(IBL) pressure to determine whether or not the clouds are 
being compressed. Thompson et al. (2004) found that two 



out of the three clouds studied have possibly been induced 
to collapse by the effect of radiation from nearby stars. The 
single system studied by |Lefloch et al.| ( |1997| was also found 
to be in this state, where the IBL pressure was greater than 



the cloud support pressure. Urquhart et al. (2006 ) performed 
a similar pressure comparison on four BRCs from the SFO 
catalogue (Sugi tani et al.|1991||Sugitani &; Ogura 1994) us- 



ing ^CO, i6 CO and C i8 (J =l->0) transitions, finding 
probable triggering in two of the clouds. Other examples are 
Morgan et ah] ( |2009| ) and |Urquhart et al.| ( |2Q09| ) where CO 
observations and signatures of photoionization were used 
to refine the northern and southern hemisphere SFO cat- 
alogues respectively, retaining only those clouds in which 



triggered star formation seems likely. Morgan et al. (2009) 



and Urquhart et al. (2009) both found that BRCs hosting 



sites of probable star formation typically had warmer exter- 
nal layers of neutral gas, approximately 20 — 30 K, compared 
to the central cloud which is at about 10 — 20 K. Urquhart 
et al. (2009| retained clouds for which a photon dominated 



region (PDR) was clearly visible, a feature which suggests 
that photoionization is taking place. These gas studies have 
been reasonably successful in identifying possible sites of 
triggering, however none has been conclusive. Additional 
evidence and tests of the accuracy of the diagnostic tech- 
niques are still required before strong conclusions can be 
drawn about the prevalence of RDI. 

A number of features have been identified in molec- 
ular line profiles that are believed to be characteristic of 
specific kinematic processes. For a gas undergoing Maxwell- 
Boltzmann thermal motions the line profile is described by 
a Gaussian distribution due to thermal broadening. Devia- 
tions from this form can give insight into the bulk motions 
of the molecular gas. A signature that is commonly inter- 
preted as representing infall comes from observations of op- 
tically thick lines such as 12 CO. If the optically thick line 
is sufficiently self-absorbed there will be two peaks in the 
line profile. For an infalling cloud, the red line profile peak 
is due to material moving away from the observer in the ex- 
terior regions of the cloud and the blue line profile peak is 
due to material moving towards the observer in the central 
regions of the cloud. Given that the central regions are at 
higher density and are therefore more likely to exceed the 
critical density for the molecular species, a blue-asymmetry 
(when the blue peak is stronger) is the expected signature 



of infall (e.g. Lee et al. 2004 Tsamis et al. 2008 Roberts 



et al.||2010| |Stahler fe Yen||2010| ). A second, optically thin, 
line such as C O is checked for a single peak to ensure that 
the two peaks from the optically thick line are from the same 
cloud rather than a superposition of two distinct objects at 
different systematic velocities. 

Although this blue-asymmetry has been observed for 
protostars (e.g. Mardones et al.|fl997 ) and pre-stellar cores 



e.g. Lee et al. 2001[ 2004) it is generally not observed in 



BRCs. Rather there is usually no clear asymmetry and some- 



times even a dominant red-asymmetry ( Thompson Sz White 
120041). For example |De Vries et al.l (120021) use the Five Col 



lege Radio Astronomy Observatory (FCRAO) to perform 
a number of molecular line observations of BRCs from the 
SFO catalogue. They found that a strong blue- asymmetry 
feature was only observed in one out of seven of the BRCs 



that they studied. |De Vries et al. (2002) proposed that this 
might be due to the shock heating the cloud from the out- 
side in, which could render the standard infall signatures 
unobservable. 

It is currently unclear what is responsible for the lack 
of infall signature in the self-absorption peaks of optically 
thick line spectra of BRCs. There are a number of proposed 
causes, for example rotation (e.g. Redman et al.|2004 ), pul- 



sation (e.g.|Keto et al.|2006|, turbulence in the core (e.g. Lee 



fc Kim 2009 Smith et al. 2012), shock heating (De Vries 



et al. 2002) or the envelope expansion with core collapse 
(EECC) model (e.g.[Keto et al |2006||Gao fe Lou|2010| [Loll 



fe Gao||20TT) |Fu et al.||2011| ). ^Thompson fe Whitej ( |2004 ) 



studied the red-asymmetric BRC SFO 11NE in IC1848 and 
also attempt to model its line profile by calculating a syn- 
thetic profile for a number of possible cloud configurations. 
They found that an EECC model gave good agreement, sug- 
gesting that SFO 11NE is in the expansion phase of RDI 



identified by Leflo ch & Lazareff| ( [1994 ). 

Understanding the reliability of molecular line diagnos- 
tics and the reason behind the lack of blue-asymmetry is 
essential if a more comprehensive picture of the effect of 
feedback and triggered star formation is to be realised. In 
|Haworth et al.| ( |2012| we tested other diagnostics of BRCs 
that use continuum and atomic line data. In this paper we 
extend this form of analysis to molecular lines: generating 
synthetic data cubes and performing standard diagnostics to 
test their accuracy and applicability, and to address sources 
of ambiguity when using them to infer whether or not trig- 
gered star formation is occurring. 



2 NUMERICAL METHOD 

We use the grid-based radiation transport and hydrodynam- 
ics code TORUS to perform non-LTE molecular line transfer 
calculations (e.g.|Harri es 2000 ; Run dle et al.|2010 Acreman 



|et al.||2010[ |Haworth fe Harries||2012| ). The details of the 



molecular line transfer algorithm are given in Ru ndle et al.| 
(|2010|). We perform a non-LTE statistical equilibrium calcu- 
lation to determine the level populations and use the result 
to calculate synthetic observations in the form of spectral 
datacubes. 

We use an accelerated Monte Carlo method (|Hogerhei-| 



jde & van der Takl 12000} to calculate the mean intensity 



in each cell. A cell-centric long-characteristic ray tracing 
scheme is used in which a number of randomly directed rays 
are traced from random locations in each cell. The frequen- 
cies are also randomly selected from a uniform distribution 
of width 4.3 turbulent line widths, centred on the rest fre- 
quency of a given molecular transition. The specific intensity 
at the end point of the ray in the cell is determined by in- 
tegrating the equation of radiative transfer along the path 
traced by the ray to the edge of the grid. The boundary 
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condition for most rays is the cosmic microwave background 
(CMB). However for the calculations in this paper there is 
a nearby O star, the effect of which has to be included. A 
ray is therefore forced from each cell towards the star, using 
the stellar effective temperature as the boundary condition 
and weighting that ray's contribution based on the assumed 
probability of it having intersected the star at random. This 
probability is simply the solid angle subtended by the star 
divided by 4ty. 

We include dust in these calculations and assume a 
canonical value for the dust to gas mass ratio of 1 x 10 -2 
in all cells that are below a temperature of 1500 K. For cells 
hotter than this we incorporate sublimation effects by set- 
ting the dust abundance to a negligible value. We assume 
spherical silicate dust grains that follow a standard inter- 
stellar medium size distribution (Mathis et al. 1977). The 
optical constants are taken fromlDraine feLee|fll984|T Given 
that BRCs are relatively young (the RHD models simulated 
200 kyr of evolution), the dust size distribution and chem- 
istry are not expected to depart much from this canonical 
interstellar medium model. This is the same dust treatment 



used in Haworth et al. 



(2012). 



Once a set of rays has been traced, the mean intensity J u 
in each cell is calculated by averaging the specific intensity 
from the rays, weighted by the line profile function 



^turb^OV 71 " 



-Av 2 /v{ u 



(i) 



where c, z^o, ^turb, and Av are the speed of light, rest fre- 
quency of the transition, turbulent velocity and velocity re- 
quired to Doppler shift z/ to v respectively. Here t>t U rb is 
imposed as 0.2kms _1 , similar to that used in Rundle et al 



( 20101 and featured in |Offner et al.| (|2QQ8|). The ra diation 
hydrodynamic models of |Haworth fc Harries (2012) (which 



are the basis for the statistical equilibrium calculations in 



this paper and are described more in section 2.1) exhibited 



strong systematic bulk motions which dominate turbulence. 
J v in each cell comprises two components, a first that is 
fixed for the cell during one set of level population itera- 
tions which describes the contribution from space external 
to the cell and a second that varies with the level popula- 
tions (which affect the source function, S u ) internal to the 
cell 



£/je- Ti ^ E^(l-e" Ti )^ 



+ 



(2) 



where V v and Ti are the intensity and optical depth along the 
i th ray. Equation [2] is solved iteratively in conjunction with 
the equations of statistical equilibrium, which determine the 
level populations and modify the source function within the 
cell 



ni 



k<l 



k^l 



^2n k A k i + (BkiJu + C k i) 

k^i 



(3) 



k>l 



where ni, Ai k , Bik and Cik are the relative fractional 
level population of level /, Einstein A (spontaneous absorp- 
tion/emission) and B (stimulated absorption/emission) co- 
efficients and the collisional rate coefficient for levels / and k 



at a given temperature. The coefficients are taken from the 
LAMDA database ( |Schoier et al.|2005| ). Initially the J = 0, 1 
relative fractional level populations are set to 0.5 and the 
other levels to 1 x 10~ 10 . 

The ray tracing and level population calculations are 
performed iteratively. Convergence is checked by comparing 
level populations from the latest and previous iteration, be- 
ing achieved where the maximum root mean square frac- 
tional difference in all levels is less than a user-specified 
value, here taken to be 1 x 10 -2 . This results in an aver- 
age fractional difference of order 10 -4 -10 -5 in the J = 2, 
1 levels which are those required for the transitions used 
in this paper. Some repeat calculations were performed to 
check that our convergence criterion was sufficient. 

A two-stage calculation is performed in which an initial 
set of iterations using rays with fixed position, frequency and 
direction is run until the level populations converge. This is 
followed by iterations using rays with random position, fre- 
quency and direction that double in number until the level 
populations again converge. The first stage of the calculation 
converges quickly, but poorly samples both the frequency 
range and the spatial extent of the grid. The second stage 
reduces the systematic and random errors associated with 
using fixed rays. This combination reduces the calculation 
time compared to using solely random rays. A typical calcu- 
lation requires around 10 iterations using a starting number 
of between 400-700 fixed rays per cell (that do not double in 
number between iterations) followed by 3-4 iterations using 
random rays which double in number with each iteration. 
We also make use of the convergence-acceleration scheme of 
|Ng| ( |1974| , which estimates an updated set of relative frac- 
tional level populations by extrapolation based on the level 
populations from the previous 4 iterations. This convergence 
acceleration is employed every five iterations. 

In this work we investigate molecular line diagnostics 
and kinematic signatures of the neutral component of BRCs. 
Due to computational expense it currently not possible to 
perform 3D radiation hydrodynamic models with chemi- 
cal evolution (Gl over et al.|[2010| . We therefore use stan- 
dard values for the molecular abundance relative to hydro- 
gen and neglect PDR and low-temperature chemistry other 
than to adopt the following conditions. An abundance drop- 
model is employed at low temperatures for CO and its iso- 
topologues to accommodate freeze-out of molecules on to 
dust grains ( | j0rgensen||2OO4| . Under this scheme molecular 
species in cells at less than 30 K and molecular hydrogen 
density greater than 3 x 10 4 cm -3 have their abundance re- 
duced by a factor of 10. The molecular abundance is set to a 
negligible value where the neutral atomic hydrogen fraction 
is lower than the conservative value of 0.1 as the gas is ion- 
ized and molecules would be dissociated (photodissociation 
codes assume no ionized atomic hydrogen, e.g. |Bisbas et al.| 
2012||Heiner fe Vazquez-Semadeni|2013| ). 



Data cubes comprising two-dimensional spatial data 
and a series of velocity channels are generated by ray trac- 
ing in a similar manner to the primary level-population 
solver, only the rays are directed towards a pixel array 
that represents the image plane at the observer position. 
The cubes are produced in units of spectral radiance, 
ergs -1 cm -2 Hz -1 sr _1 , usually simply referred to as the 
monochromatic specific intensity I u and converted into a 
brightness temperature Tb using the Rayleigh- Jeans approx- 



4 T. J. Haworth et al 



imation 



Pseudocolor] 



Var: logRho 
-19.00 - 



where c, v and are the speed of light, frequency of obser- 
vation and Boltzmann constant respectively. Extensive test- 
ing of the molecular line transfer calculations is included in 
Rundle et al.| ([2010). 



2.1 The radiation hydro dynamic models 

The density, temperature and velocity distribution that pro- 
vide the basis for the statistical equilibrium and simulated 
observation calculations are taken from the final grid states 
of the radiation hydrodynamic RDI models of |Haworth fc| 



Harries (2012) that were also used in Haworth et al. (2012) 



In Haworth &; Harries] ( |2012| , the models initially con- 



sisted of a Bonnor-Ebert sphere (BES) at the centre of the 
grid with a plane parallel ionizing radiation field impinging 
upon the left hand edge of the grid. We considered three 
different distances of the BES from the star responsible for 
the plane parallel radiation field which, due to the varying 
levels of flux incident at the left hand edge of the grid, we 
labelled the 'low', 'medium' and 'high' flux models. The ra- 
diation field parameters are all given in Table [TJ During the 
radiation hydrodynamic calculation, an ionization front was 
established which drove into the BES, accumulating a dense 
shell of material and changing the BES morphology. The 
outer layer of the shell was also photoionized and ejected 
in a photo-evaporative outflow. The manner in which com- 
pression of the BES proceeded was found to be dependent on 
the distance of the star, in agreement with previous models 



such as |Gritschneder et al. (2009) and Bisbas et al. (2011). 
We also found that inclusion of diffuse field radiation could 
significantly modify the result of the calculation. 

In this paper the final states of the most sophisticated 
RDI models, those which included the diffuse field in |Ha-| 
worth & Harries| (2012) and have been subject to a full 



photoionization and thermal balance calculation in Haworth 
et al.| ( [2012| , are used. A slice through the logarithmically 



scaled density distribution for each model at 200 kyr (the 
simulation end time of the RDI models) is given in Figure 
[I] Also included are velocity vectors and a contour corre- 
sponding to the point at which the neutral atomic hydrogen 
fraction is equal to 0.1 (above which molecular species are 
able to survive, see section [5]). Where the contour does not 
trace the dense gas, for example in the wings of the high flux 
model (the bottom panel of Figure [T]) the additional cooling 
from the full photoionization and thermal balance calcula- 
tion in Haworth et al.| ( |2012| has moved the ionization front. 
This relocation of the ionization front has a negligible effect 
on the simulated molecular line diagnostics in this paper 
since it is in regions away from the main cloud (the object 
of study) that are at low density and therefore low intensity 
relative to the cloud. 



3 SIMULATED OBSERVATIONS 
3.1 Choice of molecular transitions 
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Figure 1. Slices through the final states of the RDI mod- 
els from Haworth and Harries (2012). The panels are from the 
'low', 'medium' and 'high' flux models from top to bottom. The 
greyscale distribution is the logarithmic density, the contour is 
that at which the neutral atomic hydrogen fraction is equal to 
0.1 and the vectors represent the velocity field. Major ticks are 
separated by 0.65 pc. 



Isotopologues of CO are among the most commonly used 
species in molecular line observations, in particular for ob- 
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Table 1. A list of key parameters from the RDI models of Ha- 
worth &; Harries (2012). 



Variable (unit) 


Value 


Description 


Rc (pc) 


1.6 


Cutoff radius of initial BES 


rimax (cm -3 ) 


1000 


Peak initial BES number density 


$l ow (cm -2 ) 


9.0 x 10 8 


Low ionizing flux 


Aow (PC) 


(-10.679,0,0) 


Source position (low flux) 


$med (cm" 2 ) 


4.5 x 10 9 


Intermediate ionizing flux 


^med (PC) 


(-4.782,0,0) 


Source position (medium flux) 


$hi g h (cm" 2 ) 


9.0 x 10 9 


High ionizing flux 


^high (pc) 


(-3.377,0,0) 


Source position (high flux) 


T(K) 


40000 


Source effective temperature 


fl(Ro) 


10 


Source radius 


L(pc 3 ) 


4.87 3 


Grid size 



serva tions of BRCs (e.g. [Lefloch et al.||1997 De Vries et al 
2002||Thompson et al.|2004||Morgan et al.|2009| ). This is be- 
cause they have a relatively high abundance and low critical 
density, making them easier to observe. We therefore choose 
to generate data cubes of 12 CO, 13 CO and C 18 0. Analysis of 
molecular line data requires the use of probes which are sen- 
sitive to different conditions in the cloud. 12 CO (J = 2—^1) is 
a line which can be optically thick, with 13 CO (J = 2—tl) and 
C 18 (J = 2—^1) being optically thinner variants. These are 
combined to trace and infer the properties of the molecular 
gas. 

3.2 Simulated instruments 

We smooth the data cubes that we calculate to a Gaussian 



beam using ACONVOLVE from CIAO v4.1 (Fruscione et al. 



2006) to a size appropriate to the half power beamwidth 
(HPBW) of the simulated instrument. We choose the beam 
size given by the Rayleigh criterion, as is the case for the 
JCMT which has a 15 m dish, resulting in beam sizes of 
around 22". For comparison, the beam size of the 12 CO 
( J — 1 — >0) transition using the FCRAO in |De Vries et al. 



(2002) is 46". Factors such as instrument and atmospheric 



noise are not included. 



4 CALCULATING THE MOLECULAR CLOUD 
CONDITIONS 

The optical depth, excitation temperature and column den- 
sity of the optically thin C 18 line can be determined fol- 
lowing the method described by Myers et al. ( |1983| and 
used by, for example, |Urquhart et al. ( 2004 ) and Morgan 



et al.| ([2009]). The source- averaged optical depth of C^O is 



determined using 



lis 

Tig 



(5) 



where T13 and Tis are the peak brightness temperatures of 
the source- averaged line profiles of 13 CO and C 18 respec- 
tively, with the background signal subtracted. In this work, 
the background signal is determined from the average value 
in the ambient Hll region. Equation [5] is solved by assuming 
that the two source- averaged optical depths are related by 
their relative abundances 



Tl3 — X 13 / 18 T 18 



(6) 



where X 13 / 18 is the ratio of 13 CO to C 18 abundances. 
This abundance ratio is usually estimated based on Galac- 
tic abundance distributions and the location of the target 
in the Galaxy (iLanger & Penzias 19901) , for example being 



taken as 10 in |Urquhart et aL ( 2006| ). Here we use the ra- 
tio of prescribed abundances from Table [2] giving a ratio of 
approximately 16 for X 13 / 18 and 30 for X 12 /i 3 . We find the 
remaining single unknown optical depth numerically, using 
a decimal search. 

The gas excitation temperature is estim ated for 12 CO 
and C 18 via the same approach used in Mor gan et al. 
(2009). The equation of radiative transfer written in terms 
of optical depth r, integrated along a path length s and with 
the background subtracted gives the intensity as a function 
of the Planck function B u , the background intensity Io and 
the optical depth 



I v (s) = (B v -Io){l-e- T ). 
This can be re- written in terms of temperatures as 
1 1 



hv 
kB 



ohv/k B T e i 



ahu/k^T cmh 



(?) 



(8) 



where Tb, T ex and T cm b are the brightness, excitation and 
CMB temperatures respectively and v is the frequency of 
radiation emitted following the molecular transition ( |Rohlfs| 
fe Wilson|1996) . Equation [8] can be rearranged for the exci- 
tation temperature to 



= Tvl In 



l + T v 



T B + T,(1 



e~ T ) e J cmb 



(9) 



where T v is set to his/k-B- If the 12 CO emission is optically 
thick, which is typically expected to be the case, the term 
1 — e~ T tends to 1 and the excitation temperature T ex of 
12 CO can be estimated using 



11.06 On 



1 + 11.06 



1 



T B + 0.192 



(10) 



The criterion that 12 CO be optically thick can be checked 
by comparing the relative intensities of 12 CO and 13 CO to 
their relative abundances. If the 12 COto 13 CO line intensity 
ratio is much smaller than the assumed abundance ratio then 
12 CO is expected to be optically thick. 

Assuming local thermodynamic equilibrium (LTE) and 
that a single temperature T applies to the whole cloud, the 
kine matic temperature is simply this excitation tempera- 
ture. Morgan et al. ( |2009| ) derived the C 18 excitation tem- 
perature of clouds in addition to 12 CO and found significant 
differences. Since 12 CO and C 18 probe different parts of 
the cloud it is not surprising that they will have different 
excitation temperatures. Typically the interior parts of the 
cloud, probed by C 18 0, are expected to be cooler. With its 
optical depth known, the C 18 excitation temperature can 
be calculated independently using equation [9j giving 



T ex = 10.54 <^ In 



1 + 



10.54(1 



T B +0.221(1 -e- r is) 



(11) 



The total column density for a given molecular species is cal- 
culated following the method given in IScoville et all (119861), 



whereby the optical depth is integrated over the line profile. 
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For CO molecules, assuming a rigid rotor and that all lev- 
els are represented by a single excitation temperature, the 
column density over all levels is 



10 



15 



20 



25 



30 



35 



N ■ 



3/cb 



/lBJi(J I +l)/fc B Te: 



T ex + hB/3k 



87T 3 B/i 2 



(JZ+ 1) (1 - e-W^Te: 



dv 
(12) 

where B and \i are the rotational constant and permanent 
dipole moment of the molecule respectively. J/ is the lower of 
the two rotational levels for the transition being considered. 
The rotational constant and permanent dipole moment of 
C 18 are 54.891 GHz and 0.11 Debye respectively. When 
using equation 1 12| in this paper we use the average value of 
Ti8 calculated via Equations [5] and [6] and use the FWHM 
of the line to replace the velocity integral. With the C 18 
column density known, the H2 column density N (H2) can 
then be found using an assumed (and, since in this case it 
is prescribed, correct) abundance of C 18 relative to H2, 
namely 1.7 x 10 -7 (the prescribed value is given in Tabled 
and was taken from |Goldsmith et al.||1997| . 

Following Urquhart et al. ( 2006| (assuming spherical 
symmetry) we estimate the cloud average number density 
71h 2 using the column density via 



tv 3 N(R 2 ) 
SR 



(13) 



where R is the cloud radius. Equation |13| is derived by inte- 
grating the column density over the assumed uniform den- 
sity sphere and dividing by the circular surface presented to 
the observer. Finally the mass of the cloud can be estimated 
using 



M c] 



(14) 



where R, mu and \i are the cloud radius, atomic hydro- 
gen mass and mean molecular weight respectively. We fol- 



low Urquhart et al. ( 2006| and use a value of \i — 2.3, which 
assumes 25 per cent abundance of helium by mass. 

In addition to the above mass calculation which assumes 
spherical symmetry and uniform density, the mass can also 
be calculated by integrating the inferred column density over 
a given solid angle. That is, the total mass of molecular 
hydrogen of over all pixels % in a given region is 



(15) 



where X(j) is the abundance of species j relative to molec- 
ular hydrogen and a\ is the pixel area. 



5 RESULTS AND DISCUSSION 
5.1 Synthetic data cubes 

5.1.1 Cube construction 

We generated 12 CO, 13 CO and C 18 data cubes from the 
results of statistical equilibrium calculations in the manner 
described in sections [2] and [3] A total of nine statistical equi- 
librium calculations were run, one for each isotopologue con- 
sidered at the three distances of the ionizing star from the 
cloud (see section 2.1 for details of the model). The maxi- 
mum velocity magnitude in the data cubes is lOkms -1 and 
each velocity channel spans 0.1 km s -1 . A summary of the 
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Figure 2. Low flux model data cube 12 CO images (greyscale) 
smoothed to a Gaussian beam representative of the JCMT and 
integrated over velocity channels. The image temperature scale 
is given by the greyscale bar in integrated brightness tempera- 
ture units (Kkms -1 ). The contours are 13 CO (top) and C ls O 
(bottom). These images are all 4.87 pc to a side. 



statistical equilibrium and datacube calculation parameters 
is given in Table [2] The abundance of each species in Table 
[2] is a constant value relative to H2 , determined through a 
literature search. 



5.1.2 Edge-on morphology 

12 CO images of the clouds for an observer edge-on to the 
BRC, convolved to the appropriate Gaussian beam size, are 
given in Figures [2] [3] and [4] for the low, medium and high flux 
models respectively. Overlaid are 13 CO and C 18 intensity 
contours. These images and contours are constructed by in- 
tegrating the data cubes using the starlink software GAIA. 
The contours are chosen to give the best representation of 
the distribution of emission throughout the BRCs. The ion- 
izing star is located off the left hand edge of the images. 

The low flux 12 CO images (Figure [2]) are dominated by 
the bright bow, with weaker emission in the wings of the 
cloud. There is also some weaker emission behind the bow 
in the cloud. The integrated 13 CO and C 18 contours trace 
the 12 CO morphology well. The peak C 18 contours extend 
further into the neutral gas away from the bright rim due to 
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Table 2. The parameters used in statistical equilibrium and datacube generation calculations 



Parameter (Unit) 


Value 


Description 


^turb (kms -1 ) 


0.2 


Turbulent velocity 




200 


Number of velocity channels 


dv (kms -1 ) 


0.1 


Span of each velocity channel 


N ■ 

1 v pix 


401 2 


Number of pixels per channel 


en 


2.5 


Angular width per pixel 


tolerance 


1 x 10- 2 


Statistical equilibrium convergence checking tolerance 


n c i8 /^H 2 


1.7 x 10 -7 


C ls O abundance (Goldsmith et al.|l997|> 


ni3 CO /™H 2 


2.7 x 10- 6 


13 CO abundance rtPineda et al. 2008 and references therein) 


ni2 CO /^H 2 


8.0 x 10 -5 


12 CO abundance (JJVlagnani et al.|l988| and references therein) 
Distance of observer 


D (pc) 


1000 



4 6 8 10 12 14 16 18 20 4 4.5 5 5.5 6 6.5 7 7.5 




Figure 3. Medium flux model data cube 12 CO images (greyscale) 
smoothed to a Gaussian beam representative of the JCMT and 
integrated over velocity channels. The image temperature scale 
is given by the greyscale bar in integrated brightness tempera- 
ture units (Kkms -1 ). The contours are 13 CO (top) and C ls O 
(bottom). These images are all 4.87 pc to a side. 



Figure 4. High flux model data cube 12 CO images (greyscale) 
smoothed to a Gaussian beam representative of the JCMT and 
integrated over velocity channels. The image temperature scale 
is given by the greyscale bar in integrated brightness tempera- 
ture units (Kkms -1 ). The contours are 13 CO (top) and C ls O 
(bottom). These images are all 4.87 pc to a side. 



the lower critical density of the line. In the radiation hydro- 
dynamic calculations instabilities arose resulting in fingers 



with dense tips in the wings of the BRC (Haworth & Har- 
ries|2 012). Those dense tips are not individually resolved in 



these data cubes due to the beam size used, rather they ap- 
pear to contribute to the more widespread emission in the 
wings of the BRC. 



The medium flux 12 CO images (Figure |U show a dense 
core at the tip of the cloud and a tail that gives the ob- 
ject a comet ary appearance. There is also fairly widespread 
emission about the comet ary object which is due to fore- 
ground material, rather than material coinciding with the 
BRC itself. A channel map of the medium flux model over 
the velocity range —7.02 to +0.05 kms -1 is given in Figure 
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Figure 5. A pseudo colour channel map of the medium flux model in 12 CO. The channels, from left to right, are at —7.02, —4.70, —2.27 
and 0.05 km s — 1 . The colourbar represents the brightness temperature scale in Kelvin. Each channel is 4.87 pc to a side. 



[5] We only include a channel map into the negative veloc- 
ity range because there is no significant visual difference in 
the corresponding positive velocity channels. At low veloci- 
ties (the right hand panel of Figure [5| most of the emission 
comes from the undisturbed material in the inner core at 
the tip of the cloud, as well as from the layers of the shell 
driving into the cloud perpendicular to the line of sight. At 
higher velocities the components of the shell driving into the 
shell along the line of sight dominate. At high negative ve- 
locities the shell from the near side of the cloud is observed 
and at high positive velocities the shell from the far side of 
the cloud is observed. 

The high flux 12 CO images consist of a BRC which has 
a fairly dim bow compared to the low and medium flux mod- 
els. The weaker extended emission behind the main cloud is 
therefore more easily visible due to the reduced contrast. 
There is also a lot of visible foreground material that is not 
directly associated with the BRC, towards the right of the 
image. The relative dimness of the high flux model com- 
pared to the other two models is because less material was 
accumulated during the high flux radiatio n hydrodynamic 



calculation in Haworth & Harries ([20121. The 13 CO and 
C 18 contours have fairly similar morphology, again tracing 
the 12 CO extent of the gas well. 

In general, although the optically thin and thick lines 
may be focused on slightly different components of the cloud 
(depending on the density structure) each line traces a sim- 
ilar extent of the cloud for all models. It is therefore only 
necessary to use the combination of lines considered in this 
section to determine the average cloud conditions. 



5.2 Edge-on line profiles 

We split the edge-on image of each cloud into a series of 
20 equally sized boxes over which we calculate the average 
line profiles. These boxes are 50 by 50, 20 by 40 and 40 by 
50 pixels for the low, medium and high flux models respec- 
tively. The box sizes are chosen to provide optimal coverage 
of the BRC. Signatures in profiles such as these are used to 
infer the kinematic behaviour of the gas by observers but 
the cause of these signatures is not always clear. The inter- 
pretation of line profile features is usually based on simple 
theoretical models. Since we have directly modelled the RDI 
process and know the thermal and kinematic conditions we 



can attempt to clarify the origin of some of these signatures. 
The 20 line profiles are shown for each model cloud across 
all three considered molecular species in Figure [5] 



5.2.1 General features of the edge-on line profiles 

At this viewing angle the line profiles typically consist of 
multiple components. The primary component is a peak of 
small width centred on Okms -1 which is due to turbulently 
broadened emission from the stationary gas both within 
the cloud and from the back/foreground material. The sec- 
ondary components are due to the swept up shell of material 
driving into the cloud, in agreement with findings from LOS 
velocity profiles calculated in Tremblin et al. (20121). 



The profiles are all symmetric about the horizontal mid- 
plane of the BRC. There is no helical structure or apparent 
rotation of the cloud as is observed in some elephant trunks 



(Carlqvist et al. 1998 2003 Gahm et al. 2006). This re- 



sult is unsurprising since the starting conditions of our RDI 
models were axisymmetric. Rotating elephant trunks form 
via instability (Schneps et al. 19801 iChauhan et al. 2011) 



or the exposure of a turbulent medium to ionizing radia- 



tion (jGritschneder et al.||2009| |2010| |Ercolano et al.||2012 
Tremblin et al.|2012| ) rather than RDI. Unless the formation 
of trunks from a collection of initial inhomogeneities gives 



rise to a different velocity structure (Ma ckey &; Lim|2010 ) it 
seems that the product of RDI of larger scale existing objects 
is the BRC, which is a distinct object from the narrower, 
relatively rapidly rotating elephant trunks. The velocities in 
most profiles are in a similar range to those identified in LOS 



velocity profiles of the RDI models of Gritschneder et al 



(2010). Some of the cloud- averaged profiles from Urquhart 
et al.| (2009) also show secondary features that span the 



velocity range illustrated here, for example SF081, which 
shows two secondary peaks separated by around 9kms _1 . 
The range in velocities between the peaks in the line profiles 



observed in Urquhart et al. (2009) are too large to be the 



result of self absorption. In our the models this large velocity 
range is due to the systematic bulk shell motions relative to 
the low velocity gas encompassed by the shell. 

Features in the line profiles due to the shell have their 
intensity determined by the density at which the line be- 
comes optically thick. The peak velocity of the shell fea- 
ture depends on the viewing angle of the observer and the 
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Figure 6. Line profiles for regions of the low, medium and high flux models from left to right columns respectively. The top row shows 
the region locations. The profiles are, from top to bottom, 12 CO, 13 CO and C ls O. 
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propagation direction of the shell. For example if the shell 
is travelling perpendicularly to the observer viewing angle 
then the velocity of the shell peak will be slower than if 
the observer is viewing the shell along its propagation vec- 
tor. In the low flux case the shell layers give rise to broad 
shoulders about the central peak of relatively low intensity 
compared to the central peak. In the medium flux case the 
shell contributions to the line profile manifest themselves as 
separate peaks, since the shell is sufficiently dense and trav- 
elling sufficiently fast for the secondary peaks to be isolated 
from the turbulently broadened emission of the stationary 
cloud. In the high flux model the shell layer is thinner and 
propagating more slowly so the distinction between the un- 
compressed cloud and shell is not as clear. The peaks due to 
the shell in these models increase and decrease in strength as 
the observer moves to different viewing angles, as discussed 
further in section [531 

The position of the boxes with the most intense line 
profile peak in the low and high flux models in Figure [6] 
do not correlate between 12 CO and C 18 0. In the medium 
flux model, the position of the most intense 12 CO and C 18 
profile peaks do correlate, being situated towards the tip 
of the comet ary object in the left most column of the mid- 
dle row. This implies that in the high and low flux model 
the optically thick and thin lines are predominantly probing 
different parts of the clouds. In the medium flux case the 
optically thick and thin line both probe the same part of 
the cloud. 



Vector 
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Figure 7. A slice through the logarithmic density distribution 
of the medium flux model. Overlaid are velocity vectors and a 
contour corresponding to a neutral atomic hydrogen fraction of 
0.1. Material external to this contour around the cloud will not 
contribute to the molecular line data cubes because the gas is ion- 
ized. Note that the vortices in the hot gas are Kelvin-Helmholtz 
instabilities resulting from shear between the driving flow and 
photo-evaporative outflow. The greyscale bar is the logarithmic 
density in log 10 (gcm -3 ) and the colour bar is for the velocity 
vectors in kms -1 . This cut is 3.9 by 2.6 pc. 



5.2.3 Envelope expansion with core collapse 



5.2.2 Edge- on asymmetries 

Another interesting feature is that asymmetries in the 12 CO 
line profiles at this edge-on viewing angle are predominantly 
red. That is, non-Gaussian features with v > Okms -1 are 
stronger than those with v < Okms -1 . An example of red 
asymmetry from Figure [6] is the central row of the medium 
flux model in 12 CO where the shells are most directly prop- 
agating towards and away from the observer. The reason 
for this red-asymmetry is the dense shell of material that is 
driving into the cloud. The BRC observed from this view- 
ing angle is a three-component system, with a central (and 
ambient) gas cloud, a near-shell propagating away from the 
observer into the cloud (red-shifted) and a far-shell prop- 
agating towards the observer into the cloud (blue-shifted). 
The optically thick emission from each component will only 
be from the closest layer to the observer at a given velocity. 
The interior edge of the driving shell is typically at a lower 
density than the exterior and central regions of the shell so 
the blue peak in the line profile is from lower density inte- 
rior gas and is therefore weaker. At the frequency of the CO 
molecular transitions considered here, dust absorption plays 
a negligible role in attenuating the observed intensity from 
the far shell layer. Asymmetries become more pronounced at 
different viewing angles. For example, as the observer moves 
to view the BRC face on the shell will be denser and moving 
more directly along the observer's line of sight and the pro- 
file will be more red-asymmetric. This is discussed further 
in section [531 



An alternative explanation for the red-asymmetry in BRC 
line profiles is the envelope expansion with core collapse 
(EECC) model, in which the cloud is a two component sys- 
tem with a collapsing core and an expanding outer shell, 
(e.g.jKeto et al.|2006l|Gao fe Lou|2010l|Lou fe Gao|2011||Fu] 



et al.|20lT ). Figure 1 7\ shows a slice through the medium flux 
model density distribution and has a contour corresponding 
to the neutral atomic hydrogen fraction being equal to 0.1 
overlaid. This figure clearly illustrates that the gas outflow- 
ing towards the observer is all too ionized to harbour molec- 
ular gas (see section The contributor to the line profile 
must be the neutral part of the driving shell, the gas con- 
tained within it and any neutral foreground material. The 

Lou & Gaol (|201 lft and 



EECC models of Gao & Lou 



(2010), 



Fu et al. (2011) describe well smaller isolated starless cores 



that are not being driven by the surroundings but following 
our result probably do not extend to BRCs and RDI, where 
the high velocity motions of the dense shell dominate the 
line profile. 



5.2.4 Comparison of line profiles with observations 



SFO 81 is a BRC studied in Urquhart et al. (2009) and has a 



triple-peaked CO profile suggesting it may be viewed edge 
on. Obtaining line profiles over smaller regions of this cloud, 
in the manner of this section, would help to confirm this. 



Urquhart et al. (2009) also present a number of other 
profiles that have similarities to those here, SFO 59, 60, 73, 
80, 81, 86 and 87 are all multi-peaked in 12 CO. Interest- 



ingly, lUrquhart et al.| ([20091 suggest that SFO 80, 81 and 
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86 are all unlikely to be triggered because a PDR is not 
readily observed at 8 /im (these observations were primarily 
made using the Midcourse Space Experiment satellite, |Price| 
et al.|2001| ). SFO 73 and 87, however, do have a visible PDR 
at 8 /im and are expected to be sites of triggering. What 
the profiles of the apparently un-triggered clouds (SFO 80, 
81 and 86) have in common, compared to SFO 87, is that 
the stronger (or only visible) of the secondary peaks is blue 
shifted (the SFO 73 line profile is too complex to compare 
with the others). If these blue shifted secondary peaks are 
due to the shell, then our results imply that the observer is 
viewing the cloud from behind with the shell moving towards 
them. As such it is less surprising that the PDR is not so 
readily visible, as it would be on the opposite side of the (po- 
tentially optically thick) cloud to the observer. SFO 87 has 
a strong secondary red peak, suggesting the shell is moving 
away from the observer and the cloud is being viewed face 
on. This is supported by the fact that the PDR is readily 
visible for SFO 87. We conclude that not viewing a substan- 
tial PDR at shorter wavelengths may not be sufficient to 
rule out triggering in a BRC. Follow up analysis with longer 
wavelength Herschel or Spitzer archival data could help to 
identify a PDR in the BRCs where one was not identified at 
8 /im. 

There are also a number of wings and shoulders iden- 
tified in the line profiles given in |Urquhart et a l. (2009), 
such as SFO 51, 55, 71 and 79 that resemble the features of 
the low flux and high flux models at this inclination and all 
models at higher inclinations. 

It should be noted that although these edge-on profiles 
best illustrate the various contributing components of the 
BRC, the form of a profile changes rapidly with viewing an- 
gle. For example, the three strong peaks of the medium flux 
profile will be dominated by a single peak due to the shell 
with non-Gaussian wings as the observer moves in front or 
behind of the BRC. That an edge-on viewing angle is com- 
paratively rare is the reason that single peaked profiles tend 



to occur more frequently in observations to date (Morgan 
et al.||2009||Urquhart et~aLl|2009 ). 



5.3 The effect of viewing angle on line profiles 

We generated data cubes for 12 CO and C 18 from —90 de- 
grees (face on to the BRC) to 90 degrees (behind the BRC) 
in intervals of 15 degrees. A schematic of these viewing an- 
gles is given in Figure [8] Due to the large volume of data 
it is impractical to replicate the overlaid grid analysis pre- 
sented for the edge on viewing angle in section [572] for each 
inclination. We therefore focus on the overall variation in 
the average line profile over the cloud. 



-90° < 




-• 90° 



Figure 8. A schematic of the viewing angle convention used. At 
—90° the observer is facing the bright rim of the class A cloud. 
At 90° the observer is behind the cloud. 



dramatically. As discussed in sectio n [5~2] this is because the 
shell is optically thin to the the C 1 O line and so the line 
profile peak comes from the interior cloud at all viewing an- 
gles. The shell is optically thick to the 12 CO line meaning 
that the line profile will change with viewing angle as the 
motion of the shell along the line of sight changes. Both the 
shell and central cloud are identified in the optically thick 
and thin line profiles, but the relative strengths differ. 

The high flux model shows similar behaviour to the low 
flux model, with the optically thick and thin peaks sepa- 
rated, but it is not as extreme. This is because the shell 
accumulated before driving into the cloud had lower mo- 
mentum than in the other models, meaning that the shell 
rapidly reached pressure equilibrium with the cloud. The 
shell only continues to propagate into the cloud due to the 
rocket motion resulting from a weak photo-evaporative out- 
flow and so the velocities are lower. 

In the medium flux model the shell is sufficiently dense 
that it is optically thick to both the 12 CO and C 18 lines. 
Therefore both line profile peaks come from a region of the 
BRC with similar kinematic properties and their peak ve- 
locities vary with viewing angle in the same way. 

Examples of 12 CO and C 18 cloud -averaged profiles at 
a viewing angle of —60 degrees are given in Figure [10] These 
illustrate the points discussed in this section, showing that 
the low and high flux models' optically thin lines stay cen- 
tred at low velocity whereas the medium flux peak moves 
to follow the optically thick line. The medium and high flux 
model optically thick and thin line peaks are only slightly 
separated at this inclination, whereas the low flux peaks are 
widely separated by 2.9 km s -1 



5.3.1 The variation of the line profile peak intensity 
velocity 

We plotted the velocity at which the average line profile over 
the cloud is at maximum intensity (hereafter referred to as 
velocity for brevity) as a function of viewing angle for both 
the 12 CO and C 18 lines across all models in Figurej9] 

In the low flux model (the top panel of Figure [9| the 
optically thin C 18 line is constant as a function of view- 
ing angle, whereas the optically thick line varies in velocity 



5.3.2 The variation of the optically thin line profile 
FWHM 

The C 18 FWHM (which is used to calculate the column 
density, c.f. equation 12) remains fairly constant at about 
0.5kms _1 in the low and high flux models. This is because 
the optically thin line profile peak is determined by the cen- 
tral cloud at all viewing angles. 

Conversely, the FWHM of the medium flux C 18 line 
exhibits a maximum at low viewing angle, decreasing by up 
to 20 per cent as the observer moves to face the object from 
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behind or face on. This is illustrated in Figure [II] The pro- 
files for which these FWHM are calculated are averaged over 
a constant number of pixels across viewing angles, centred 
on the area of peak emission and not diluted by the ambient 
medium. This variation in FWHM is hence not due to vary- 
ing the size of the region over which the profile is averaged, 
modifying the size of the line profile peak. Rather, the reason 
for this variation in the FWHM is that the shell is optically 
thick to C 18 0. As the observer moves to higher viewing an- 
gles a single, denser component of the shell than that seen 
edge-on dominates the profile. There is also a smaller distri- 
bution of velocities about the peak since only a single shell 
is contributing to the line profile rather than two. The result 
is a slightly stronger peak that has a smaller FWHM. This 
variation in the line profile is illustrated in Figure [l2| where 
the medium flux C 18 profile is shown as the observer moves 
from degrees to —90 degrees in 30 degree intervals. 



5.3.3 The variation of the line profile symmetry 
parameter 

The symmetric nature of a profile can be quantitatively ex- 
pressed using the profile symmetry parameter SV based on 
an optically thick and an optically thin line, defined as 



SV 



Vthic 



V thh 



(16) 



where Vthick, Vthin and AVthin are the source- averaged spec- 
trum peak velocities of the thick and thin lines and the 
FWHM of the optically thin line respectively ( |Mardone"s| 
et al.|1997| ). This value gives an indication of the asymmetry 
in a line profile, with negative values blue-asymmetric and 
positive values red- asymmetric. pVlardones et al. (|l997) sug- 
gest that values in the range —0.25 < SV < 0.25 should be 
considered symmetric. This is usually applied to a single op- 
tically thick line with self absorption to determine whether 
the red or blue motions are predominantly self absorbed. 
The optically thin line would have a similar linewidth as the 
thick line and typically \SV\ < 1. In this paper, although 
both the optically thin and thick lines come from the BRC, 
they sometimes probe different regions (the shell and the 
cloud behind the shell) and so the optically thin and thick 
line peaks may be located at different positions giving rise to 
larger values of SV. Equation [Hp] is therefore more of a peak 
separation function than a symmetry function in this paper, 
though we still refer to profiles as symmetric or asymmetric 
depending on the value of SV. We calculate SV for all three 
clouds over each viewing angle to see if there is a systematic 
variation. The results of this analysis are given in Figure p^] 
The low flux and high flux models show a systematic 
transition in SV with viewing angle. For viewing angles 
where the observer is looking face on to the BRC (< de- 
grees using the convention in Figure |8j), SV shows that there 
is no or strong red-asymmetry. When the observer is behind 
the BRC (> degrees using the convention in Figure [8| SV 
shows that there is no or strong blue- asymmetry. As already 
discussed, this is due to the near shell motion (probed by 
the optically thick line) relative to the interior cloud motion 
(probed by the optically thin line). The low values of SV 
are similar to those found in, for example, |De Vries et al.| 
(2002). The large values arise when the shell peak becomes 



tions of SV have not yet been performed for BRCs where 
the optically thick and thin lines are widely separated and 
would therefore give rise to larger values. Some of the pro- 
files given in Urquhart et alj (|2009| visually suggest a large 



separation between the optically thick and thin line peaks, 
such as those in SFO 71 and 73. 

In the medium flux model there is no clear transition 
between the dominance of red or blue- asymmetry. This is be- 
cause the shift in velocity demonstrated in Figure[9]is similar 
for the optically thin and thick lines since they both trace 
the dense shell region of the cloud. Red and blue- asymmetric 
profiles can still arise in the medium flux model (for exam- 
ple at 15 and —45 degrees in Figure 13), though typically 



the profiles are symmetric. The values of SV obtained for 
the medium flux cloud are very similar to those obtained for 
type B-C BRCs inlDe Vries et al.l (12002). 



Given the above, the only way in which a BRC will have 
a blue-asymmetric profile is if it is viewed from the rear and 
the shell is not sufficiently dense to dominate the profiles of 
both the optically thick and thin lines. Requiring that the 
BRC be viewed from behind (with a viewing angle > de- 
grees using the convention given in Figure [8} would reduce 
the chances of observing a blue-asymmetric line profile by 
50 per cent. This is further reduced depending on the shell 
densities of real BRCs. These results do therefore provide 
an explanation for the lack, but not complete absence of, 
blue-asymmetric BRCs and rather a dominance of symmet- 
ric profiles. 

5.4 Molecular cloud conditions 

We applied the diagnostics detailed in section [4] to each 
of the clouds to calculate the mass, temperature and col- 
umn density at an edge-on inclination (0 degrees using the 
convention given in Figure [8]). The brightness temperature 
of the cloud that is used in equation [5] is the peak of the 
source- averaged spectrum with the background signal sub- 
tracted. We used GAIA to obtain a single averaged spectrum 
for each BRC and fit it with a Gaussian profile to obtain the 
peak value and FWHM. The radii of the clouds in the low, 
medium and high flux models were estimated to be 0.51, 
0.19 and 0.8 pc respectively from their spatial extent in the 
simulated images. 

The inferred properties are all presented in Table [3] 
along with the conditions from the model grid for com- 
parison. The mass calculated assuming spherical symmetry 
(equation 14) is given by M sp h and the mass calculated by 



integrating the column density (equation 15) is given by 

Mint. 

The optical depths are similar to those found obser- 



vationally by, for example, Urquhart et al. (2006), Morgan 
et al.[fl2009| ) and given in |Urquhart et al.| ( |2009| ~Other than 
for CO in the low flux model, the excitation temperatures 
are consistent underestimates of the prescribed neutral gas 
temperature of 10 K. It is often assumed that the BRC is 
in LTE and that therefore the 12 CO excitation tempera- 
ture can be used to describe the kinetic temperature of the 
cloud. The results here suggest that this could be inaccu- 
rate by up to a factor of 1.6. The column densities are also 
similar to those found observationally, for example Morgan 
et al.| (|2009t and |Urquhart et al.| fl2009|). Our values are 



stronger than the low velocity central cloud peak. Calcula- 



slightly lower than those from Urquhart et al. (2006), due 
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Table 3. The cloud conditions calculated using 12 CO, 13 CO and C ls O. The columns are, from left to right: which model the conditions 
have been calculated for, the optical depth of 13 CO and C ls O, excitation temperatures for 12 CO and C ls O, the column density calculated 
from synthetic observations, the column density from the model grid, the C ls O FWHM, the cloud mass calculated assuming spherical 
symmetry, the cloud mass calculated by integrating the column density over the cloud and the known cloud mass from the model grid. 



Model 


T13 


T18 


T12 


Tig 


log(7V(iJ 2 )) 


grid \og(N(H 2 )) 


FWHM, Ai/is 


M sph 


Mi nt 


M gr id 








(K) 


(K) 


(cm" 2 ) 


(cm" 2 ) 


(kms- 1 ) 


(M ) 


(M ) 


(M ) 


Low 


2.81 


0.18 


12.6 


6.4 


21.11 


21.16 


0.47 


100 


19 


19 


Medium 


2.56 


0.16 


6.4 


4.5 


21.61 


21.37 


1.88 


45 


9 


4 


High 


1.72 


0.11 


6.6 


6.9 


20.99 


20.82 


0.60 


185 


36 


21 



to the higher 12 CO excitation (and hence kinetic) tempera- 



tures that they obtain, of order 30 K. Urquhart et al. (2006) 



attribute this to some internal heating mechanism such as a 
young stellar object (YSO) or ultra compact (UC) Hn re- 
gion which are not present in our models. Rather, for fully 
neutral gas we prescribe a minimum temperature in the pho- 
toionization calculation of 10 K. The inferred column densi- 
ties correspond reasonably well to the column density from 
the model grid, agreeing to within 11, 43 and 32 per cent 
for the low, medium and high flux models respectively. 

When assuming spherical symmetry, all of the inferred 
masses are larger than the actual mass in the region over 
which the diagnostics were performed. The discrepancy 
ranges from a factor of 5 to 11 under this assumption. For the 
integrated column density method the agreement is much 
better, with agreement to within one solar mass in the low 
flux case up to a factor 2.25 in the medium flux case. 



In Haworth et al. (2012) we calculated the mass and 



temperature of the same clouds using greybody fitting of 
the cloud SED. Comparing to the results here to those from 



from Haworth et al. (2012), calculating the cloud tempera- 



ture based on greybody fitting of the system SED is a more 
accurate technique than the molecular line diagnostics, typ- 
ically agreeing to within 1-2 K. This is because the former 
diagnostic is based on more information, over a larger fre- 
quency range from the cloud and also makes fewer assump- 
tions in converting observational intensities to a tempera- 
ture. For the SED fitting temperature diagnostic, the main 
assumptions are an index which describes the frequency de- 
pendency of dust emissivity and that the SED can be fitted 
as a greybody. This diagnostic does also probe the whole 
cloud. However in the molecular line temperature diagnostic 
(see section [4] equations [5] through [TT|) assumptions include 
that the optical depths of 13 CO and C 18 can be related 
by their abundances (equation |6| , that the cloud is in LTE 
and that a single temperature applies to the whole cloud. A 
single line will also only give a diagnostic temperature for 
the subset of the cloud that it probes. The masses calculated 
using molecular line diagnostics are more accurate, differing 
at most by a factor of 2.25 in the integrated column density 
method compared with a difference of up to a factor of 4 
via SED fitting. This is because SED fitting assumes a con- 
stant dust to total mass conversion factor between different 
BRCs. 

A measure of stability against collapse of a BRC is 
given using the virial theorem, comparing the IBL and neu- 
tral cloud pressures (Hartmann 2009 Haworth et al.|20"l2 ). 
In general, the cloud masses have been overestimated here. 
Given this, the neutral cloud pressure and hence the stability 



against collapse may also be overestimated when using neu- 
tral cloud properties based on molecular line calculations. 



6 SUMMARY AND CONCLUSIONS 

We have generated synthetic molecular line observations 



of the models of RDI from Haworth & Harries (2012) 



1 



Using data of the CO, 13 CO and C i8 (J = 2 
transitions we have analyzed line profiles over the imaged 
BRCs and replicated standard diagnostics to calculate the 
BRC properties. Using the derived conditions and line 
profiles we have searched for signatures of RDI and tested 
the accuracy of the diagnostics. We have also investigated 
the variation of BRC line profiles with observer viewing 
angle. We draw the following main conclusions from this 
work: 

1. The synthetically imaged BRCs have a similar 
morphology to real BRCs. The optically thin and thick line 
integrated intensities all trace a similar extent of the cloud 
in each model. 

2. The lack of blue- asymmetry observed in BRC line 
profiles can be explained by the shell of material that drives 
into the cloud. If the shell is very dense then it may be 
optically thick to both 12 CO and C 18 0. If this is the case 
then the profiles of both lines are dominated by emission 
from the shell and have very similar peak velocities that 
result in a symmetric profile. In the intermediate case 
when the shell is less dense the optically thick line profile 
is dominated by the high velocity shell and the optically 
thin line dominated by the low velocity cloud interior to 
the shell. This results in an asymmetric line profile. For 
asymmetric profiles, when the observer is facing the BRC 
the shell is moving away into the cloud and there will 
be a red asymmetry. If the observer views the BRC from 
behind then the motion of the shell will be towards the 
observer and there will be a blue asymmetry. If the shell is 
sufficiently weak then it will not contribute to the profile 
and is likely that RDI will not be occurring. 

3. By examining the known motion of material in 
the neutral gas from the model grid we rule out envelope 
expansion with core collapse (EECC) as the cause of the 
asymmetry in the simulated line profiles. This is because 
expansion from the outer layers of the BRC towards the 
observer (a key feature of the EECC model) is from gas 
that is ionized, meaning no molecular gas exists in these 
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regions and they cannot contribute to the line profile. 

4. The profiles that we obtain exhibit shoulders and 
wings that resemble observations (see Figure [6]). At edge-on 
viewing angles both the near and far shell, as well as 
the gas interior to the shells, contributes to the profile. 
This gives rise to more complex profiles with up to three 
peaks. That such complex profiles exist in observations 
to date, for example the profiles of SFO 59, 60, 73, 80, 
81, 86 and 87 from Urquhart et al. (2009), is evidence 



of a shell contributing to the line profile. These systems 
should be investigated more closely using spatially resolved 
profiles. At other inclinations the profile is typically either 
invariant (for the optically thin line) or becomes dominated 
by a single peak due to the shell with non-Gaussian 
wings (for the optically thick line). Such profiles are 
most common in observations due to the higher proba- 
bility of viewing a BRC at an inclination that is not edge on. 

5. For BRCs, failing to identify a PDR at shorter 
wavelengths does not necessarily rule out RDI. If the cloud 
line profile has a secondary strong blue peak then the 
shell may be driving towards the observer, something that 
(according to the models here) only happens if the observer 
is behind the BRC. As such the PDR is on the opposite 
side of the cloud so may be more difficult to detect if the 
foreground cloud is optically thick. Examples of this could 
be SFO 80, 81 and 86 which were identified as not likely 
being sites of triggering in Urquhart et al. ( 2009 ) at 8 /im 



(primarily using Midcourse Space Experiment data), but 
have secondary blue peaks. Conversely SFO 87 (with a 
secondary red peak that suggests the shell is moving away 
from the observer and the BRC is viewed face on) does have 
a PDR identified. Analysis of the clouds where no PDR was 
detected at 8 /im using longer wavelength data such as that 
taken with Herschel or Spitzer may help to identify a PDR. 

6. The cloud conditions that we infer by replicating 



the diagnostics of, for example |Urquhart et al. ([2006 ) and 



Morgan et al. (2009), yield results that are similar to those 



found observationally. The inferred kinetic temperature 
differs from the prescribed temperature by up to a factor 
of 1.6. The column densities for low, medium and high flux 
models agree with those from the model grid to within 11, 
43 and 32 per cent respectively. The cloud masses calculated 
assuming spherical symmetry are overestimates by up to a 
factor of 11. Integrating the column density over a region to 
determine the mass yields much more accurate results, at 
worst differing from the grid mass by a factor of 2.25 and 
agreeing more closely for the other models. By comparing 



with the results from Haworth et al. (2012) we conclude 



that calculation of cloud temperatures via greybody fitting 
of the SED is more accurate. However, the mass calculation 
is more accurate using molecular line diagnostics because 
the SED fitting assumes a constant dust to total mass 
conversion factor between clouds. 



This paper marks the third in a series ( Haworth &; Har- 
ries 2012; Haworth et al.|[2012 ) that have attempted to re- 
solve some of the disparity between models and observations 
of RDI by testing assumptions in the models and simulating 
observations. At present there are still a number of untested 



approximations in radiation hydrodynamic modelling, for 
example the treatment of metals and photodissociation. Star 
forming regions also have much more complex morphology 
than the simple case of isolated RDI, which is now being 
modelled (e.g. |Dale fe Bonnell|2012| |Walch et al.|2012| ). Fu- 
ture simulated observations of these more complex geome- 
tries may reveal and resolve further difficulties in translating 
between theory and observation. We next intend to build on 
the work in this paper by using far-infrared observations of 
BRCs to try and identify triggering in the clouds that we 
predict are being viewed from behind. 
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Figure 9. The variation of the velocity at which the average 
line profile over the cloud is at maximum intensity with viewing 
angle for the low (top), medium (middle) and high (bottom) flux 
models. 
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Figure 10. 12 CO and C ls O line profiles at a viewing angle of 
—60 degrees for the low, medium and high flux models from top 
to bottom. The vertical lines run through the optically thin line 
profile peak. 
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Figure 11. The variation of the medium flux model C ls O 
FWHM with viewing angle. 
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Figure 12. The variation in C ls O line profile with viewing angle 
for the medium flux model. 
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Figure 13. The variation of the profile symmetry with viewing 
angle for the low (top), medium (middle) and high (bottom) flux 
models. The asymmetry parameter SV (equation |16| is deter- 
mined by the difference between the lines of Figure [9] divided by 
the FWHM of the C ls O line. 



